Shear viscosity and spin diffusion in a two-dimensional Fermi gas 
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We investigate the temperature dependence of the shear viscosity and spin diffusion in a two- 
dimensional Fermi gas with contact interactions, as realized in ultra-cold atomic gases. We describe 
the transport coefficients in terms of a Boltzmann equation and present a full numerical solution for 
the degenerate gas. In contrast to previous works we take the medium effects due to finite density 
fully into account. This effect reduces the viscosity to entropy ratio, rj/s, by a factor of three, and 
similarly for spin diffusion. The trap averaged viscosity agrees well with recent measurements by 
Vogt et al. [Phys. Rev. Lett. 108, 070404 (2012)]. 
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I. INTRODUCTION 

Ultra-cold atoms have emerged as a versatile system to 
study quantum effects in strongly interacting fermionic 
and bosonic many-body systems with excellent control 
over the Hamiltonian parameters Transport proper- 
ties provide particularly valuable probes which can re- 
veal the nature and strength of the effective interaction 
between particles. The shear viscosity 77, e.g., measures 
the internal friction in a quantum fluid, which is lowest 
for strongly interacting systems. For certain relativistic 
gauge theories the ratio of the shear viscosity to the en- 
tropy density s has been computed using the AdS/CFT 
correspondence and takes the value (r]/s) m i n = %j (^irks) 
2}. It has been conjectured that this value provides a 
lower bound also for a wider class of relativistic field theo- 
ries [U , and quantum fluids which saturate this bound are 
denoted as "perfect fluids" |J|. Subsequently quantum 
fluids ranging from (non-relativistic) ultra-cold atoms to 
(relativistic) quark-gluon plasmas have been investigated 
in the search for a perfect fluid which saturates this 
bound [4]. In the solid state context, the viscosity of 2d 
graphene layers has been shown to decrease logarithmi- 
cally with increasing temperature Q coming reasonably 
close to the lower limit. Another example is the vis- 
cosity of the unitary Fermi gas in three dimensions (3d) 
which has been measured recently and comes rather 
close to the hypothetical bound for temperatures below 
the Fermi temperature. This is in agreement with calcu- 
lations based on kinetic theory for low [qI-FTTI) and high 
temperatures [L2T4i5| ]. These calculations have been con- 
firmed and refined in approaches based on the Kubo for- 
mula with self-energy [l6| and full vertex corrections [TtJ 
and recently also in the form of a Quantum Monte Carlo 
simulation [18j . A similar lower bound is also seen in the 
spin diffusion coefficient D which has a minimum close to 
the quantum limit ~ K/m [H, H(| again in good agree- 



ment with calculations based on kinetic theory |19l , |21| . 

Recently, interacting ultra-cold gases have been real- 
ized in two dimensions (2d) where quantum and inter- 



action effects are even stronger than in three dimensions 
[22h]26| . Measurements for a trapped two-component 2d 
Fermi gas with strong interactions have found the viscos- 
ity to decrease with decreasing temperature and increas- 
ing interaction strength [27]. 

In this work we compute the shear viscosity r\ and 
the spin diffusion coefficient D of an interacting two- 
component 2d Fermi gas within kinetic theory. Previ- 
ous studies have investigated transport without medium 
effects on the scattering cross section j28mfj and found 
a minimum value 77/s 20(?7/s) m i n [28|, |29(. We now 
include medium scattering, which is known to strongly 
influence the dynamical properties [l3l - [l5l . I23I l3l| . and 
find that it substantially lowers the viscosity by a factor 
of about three already above T c . For the spin diffusion 
coefficient we find a similar reduction. 

The organization of the paper is as follows: In section 
HIl we introduce the model Hamiltonian, the T-matrix in 
medium, and then derive the quantum kinetic equations 
in section IHII A discussion of the zero mode in spin dif- 
fusion is found in section MI CI Readers familiar with the 
Boltzmann approach may skip ahead directly to the re- 
sults which are presented in section IIVI We close with a 
comparison to experiment in section [V] and conclude in 
section IVII 



II. THE MODEL 

We consider two species a = f, I of fermionic atoms in 
two dimensions, which are described by the grand canon- 
ical Hamiltonian 



ka 



krj 
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with the free single-particle dispersion Ska = k 2 /2m cr 
(H = 1), spin-dependent chemical potential \i a and area 
V. The model is formulated for the general case of a 
heteronuclear mixture with different values for rzif, toj,, 
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fXf, and fj,i, however all numerical calculations are car- 
ried out for the balanced case fif = ^ = /x for equal 
masses m — m-f = mj, in view of the experiment [27]. 
At ultra-cold temperatures the attractive s-wave contact 
interaction go acts only between different species due to 
the Pauli principle. The two-body scattering between a 
single f and 4- fermion is given by the exact two-body 
T-matrix [H [H 



T (E) = 



2ir/m r 



Hsb/E) 



(1) 



in terms of the reduced mass m~ = + . The 
pole at E = — eb < corresponds to the two-body bound 
state, and the binding energy eb = l/(2rn r a| D ) defines 
the 2d scattering length ■ This bound state is always 
present in an attractive 2d Fermi gas [H, [34{ . The vac- 
uum scattering amplitude for two particles with momenta 
k and — k in the center-of-mass frame is then given by [lj 
f(k = |fe|) = 2m r T {k 2 /2m r ) = 4n / [\n(l / k 2 a 2 D ) + in}. 
The scattering amplitude depends logarithmically on en- 
ergy in both the low- and the high-energy limits: this is 
due to anomalous (logarithmic) quantum corrections to 
the classically scale invariant contact interaction [35], [36[ . 

At finite density the two-particle scattering in the pres- 
ence of the medium is described by the many-body T- 
matrix T(q,u>). It can be calculated from the solution 
of the Bethe-Salpeter equation for the ladder approxima- 
tion of repeated particle-particle scattering [37H39| , and 
in the general case of spin imbalance it is given by [23| 



T 1 (q,uj)=% 1 (uj + iO + /x t + H\.-w q ) 
d 2 k /°(fc) + /°(fc + g) 



(27r) 2 U) + iO + /Xf + (JL± - Efef - Ek+ql 

with the Fermi-Dirac distribution 

f 



(2) 



(3) 



j3 = l/(fcsT) and u q — q 2 /(8m r ). While this integral is 
known analytically at T — [231 ]. at finite temperature 
we can perform only the angular average analytically but 
have to compute the radial integral numerically. Com- 
pared to the case with the bare T-matrix this increases 
the numerical effort in solving the Boltzmann equation 
considerably. 



III. TRANSPORT PROPERTIES FROM THE 
KINETIC APPROACH 

We use the kinetic approach to derive the transport 
coefficients in our system. This approach is valid pro- 



vided quantum interference effects are negligible and de- 
viations from wcll-dcfincd quasiparticles are small, which 
we assume in the following. This assumption is question- 
able for temperatures well below the Fermi temperature 
Tp and the results should be compared to calculations 
within a formalism which does not require the quasipar- 
ticle picture to be valid [13, [3 ■ 
The Boltzmann equation reads 



[d t + Vd x + F ext d k ]f a (k) = -icolll/er, f-a] 



(4) 



which is an intcgro-differential equation for the quasi- 
particle distribution function f^^(k). The left-hand side 
accounts for perturbations driving the system away from 
the equilibrium situation, while the right-hand side ac- 
counts for collisions between quasiparticles. 



A. General formalism: variational approach 

The approach we take is standard but we present it 
such that generalizations are possible in a straightfor- 
ward manner. An excellent account of this approach has 
been given in Refs. 40-42] among others. The left-hand 
side of Eq. (fj| consists of three independent differential 
operators and is henceforth referred to as the driving 
term, owing to the fact that they drive the system away 
from equilibrium. The individual terms describe tem- 
poral variations (dt), spatial variations (d x ), as well as 
external forces (dk), while the right hand side describes 
collisions due to interactions (or in other systems also 
disorder) and consequently is called the collision inte- 
gral. One can solve for the non-equilibrium distribution 
function in the linear response regime assuming that the 
deviation from the equilibrium distribution function can 
be obtained in an expansion in the perturbation. This 
schematically assumes the form 

f a (k) = f a (k) + i/°(fc) (I - f a (k)) fl(k) (5) 

for fermions where f®(k) is the Fermi-Dirac distribution, 
see Eq. ([3]) , and /* is linear in the perturbation and other- 
wise a generic function (this is true for any type of pertur- 
bation considered here). The factor ~f®(k) (l — f®(k)) 
is introduced for later convenience. In this limit it is 
consistent to approximate the collision integral by 



Icoll[fa, f-tr] 



C[f 1 a ,f 1 _ a ] + 0(f 2 ,f 2 _ a ) 
C\fl,fK\ (6) 



with 
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f-a] — 7p 



fcl,q 



5{£ka + fffci-cr — £k+qa — 



T(k + fei,£fc CT + e kl - a - fl a - fl- 



x [f° (k) f_ a (fcj) (1 - /° (fc + q)) (1 - /° CT (kx - q))] [fl (k) + f\ (fex) - (fe + q) - /i ff (fc x - g)] (7) 



where J fc = /(0. 

We are interested in the stationary solution, i.e., 
<3*/<x = on the left-hand side of Eq. Qj. To linear 
order in the perturbation one can replace f a — > f® on 
the left-hand side and write 



(8) 



where D a in the most generic case is a tensor differential 
operator acting on /° and a labels the perturbation we 
consider. In general we have 



T 



(i - /°) W 



(9) 



where we use the Einstein summation convention. At 
this point we have introduced F£? as a generalized force 
field and Pj as a generalized projection. For reasons of a 
concise presentation we assume from now on that we can 
absorb the spin-dependence of F % J into the factor Pj and 
work with F tJ only which acts in the same way on both 
spin species. For concreteness, in the case of an electrical 
conductivity we have F % J = E l 5ij and Pj — ev t ka Sij. 
This general form also dictates the form of the ansatz for 
which we choose as 



fl{k) = F i J^{k) = F i Jl i Jg a {k). 
The generalized current then reads 



(10) 



= E / 

= E^ / fio- 
= -E/x^s 

„ J k 



Ik 

-(x ij p 



(11) 



where \x lJ ) — (x?)Xi ) is a spinor and the components 
are themselves vectors in function space. In the last line 
we have introduced a scalar product. Using this defini- 
tion of a scalar product we can also define 



C[x ij ] = \{x ij \C\ X i3 )F ij . 
We can now introduce a functional 

Q[x lJ } = s[x l3 ]+c[ x ^} 

whose extremum in function space 



d x Mx l 







(12) 
(13) 
(14) 



can be shown to lead to the Boltzmann equation for the 
respective species. Conversely, the Boltzmann equation 
implies that the current reads 



fi = -S[ X f ma *}F' 11 = 2C[ x ^ max ]i^' 
= -2Q[ X f max }F l i . 



(15) 



The proper strategy to solve the Boltzmann equation is 
thus to maximize the functional Q[Xa ] l° r X % J = I l J 9a{k) 
by varying g a {k). This is done by identifying the phys- 
ically most relevant modes g n a(k) and writing g a {k) as 
an expansion with respect to these modes: 



9<r{k) = E ^n9n<r{k) 



(16) 



Maximizing Q[g a (k)] with respect to the expansion coef- 
ficients A„ leads to a matrix equation for A n which can 
be solved by matrix inversion. Usually the most relevant 
modes are the slow modes which are related to almost 
conserved quantities whose relaxation is described by the 
collision kernel. 



B. The shear viscosity within Boltzmann theory 

We consider a two-component Fermi gas in its most 
general form, allowing for different chemical potentials 
for the two species, i.e., ^ and fi^ and a species depen- 
dent mass m a . We are concerned with a system with- 
out external forces, i.e., -F xt = in its stationary state 
dtfa = 0. We assume a uniform flow in x-direction and a 
velocity gradient in y-direction, i.e., u — (u(y),0) which 
leads us to analyze the Boltzmann equation according to 



vd x f a (k) = -icoil[/ CT ,/- CT ] 



(17) 



The collision term for the contact interaction in its lin- 
earized version was introduced in Eq. ([7]). The driving 
term reads 

' m a i ay 



Following the logic of section MI Al we define more gen- 
erally 



F X J = diUj + djUi - -6 lj d l ui 



(19) 



with Ui being the components of the flow velocity of the 
fluid. The generalized current is the viscous part of the 
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stress tensor describing hydrodynamics in two spatial di- 
mensions, 



3 



(20) 



where r\ is the shear viscosity and £ the bulk viscosity. 
Combining (f20|) and (fTTj) one obtains 



(21) 



for the exact solution \x %3 )- The variational principle 
provides us with a lower bound. If we make an ansatz 
I Y ansatz ) using a finite function set g n <r{k) this implies 

m 



V > s[ x a 



(22) 



where |x^ax tz ) corresponds to the optimal choice for a 



finite number of the parameters X n introduced in Eq. (|H: 
which maximizes Eq. (|14[) . In the case of the viscosity 
there is no conserved quantity which is excited. We found 
that, just as in the three-dimensional case [16j . the choice 
for the modes Eq. (HI 



gAk) = i 



(23) 



yields results which are very close to the exact result. We 
have checked this statement for different sets of modes, 
for instance g n <y{k) — k n for n = 0, . . . , N up to N = 10 
as well as Chebyshev polynomials up to the same order 
and have found no pronounced differences. 



C. Spin diffusion within Boltzmann theory 



The spin conductivity is again bounded from below by 



> 5[^ ansatz ] 



(27) 



and we can deduce the spin diffusion coefficient D via 
D — (28) 



1± 

Xs 



with the spin susceptibility of the free Fermi gas 



Xs 



t /°(fc = 0) + mjfjk = 0) 
2tt 



(29) 



In the case of the viscosity the driving term does not 
couple to a conserved quantity such as the total energy or 
the momentum. Consequently, the variational approach 
can be employed with relatively little care and very few 
modes suffice to solve the problem essentially exactly. In 
the case of the spin diffusion this ceases to be true and 
the driving term in general does not decouple from the 
momentum mode. The momentum mode corresponds to 
the choice 



9<y 



(30) 



and if we calculate the overlap of the momentum mode 
with the driving term within this variational ansatz it 
reads 

( X \V S ) = - (rr^ ln(l + e^t) _ m;ln(l + e^)) . 



Spin diffusion in a metal describes the response of a 
system of fermions to a gradient in a magnetic field. In 
our setup this translates to the two fermion species re- 
sponding to gradients in chemical potentials, which are 
opposite for the two species. Again we discuss the most 
generic situation, which is that there are two species of 
fermions with different chemical potentials fi-f ^ /zj, and 
different atomic masses ^ m^. We assume there is a 
chemical potential for the individual atoms fi a + r ■ V/x CT . 
The distribution function is accordingly driven out of 
equilibrium by 

^ = -^^/°(l-/°) • (24) 

In the following we assume that the absolute value of the 
gradient is the same for both species but counteracts, 
V/z CT = crV/x, such that 

VI = (1 - /°) . (25) 

Again, we identify the generalized force and projector 
F ij = i ),,„),, 

I*J = ™i,Jn ■ ( 26 ) 



This is zero if /x-j- = /i| = /i and = = m, meaning 
the momentum mode is not excited. If these conditions 
do not hold the momentum mode is excited and it can- 
not be relaxed. This formally leads to an infinite spin 
conductivity a s . In metals the standard situation is spin 
balance with a finite spin conductivity, as has been dis- 
cussed recently in the context of graphene [43j . In the 
experiments under discussion two clouds of different spin 
species are prepared to collide in the center of the trap. 
If the two clouds are equal in number of particles and 
masses the unified cloud will reside in the center of the 
trap. One could excite the zero mode if one prepared dif- 
ferent densities and/or different masses for the different 
spin species. The zero mode of the spin diffusion then 
has a very simple and intuitive physical meaning and it 
corresponds to a center of mass motion. 

In our concrete setup in a balanced system we work 
with the choice 

9a = m (31) 

which is not a zero mode of the collision integral and 
has finite overlap with the driving term. We have again 
checked more generic mode choices and found this to pro- 
vide an excellent variational ansatz. 
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IV. RESULTS 

Wc have obtained the viscosity and spin diffusion from 
the variational approach using the variational ansatz 
functions introduced in Eq. ([25]) and Eq. (l3"Tj) . The trans- 
port coefficients are normalized by the respective thermo- 
dynamic quantities density, pressure, and entropy den- 
sity, and for consistency they all have to be computed at 
the same level of approximation. A definite prescription 
is provided by the large- N expansion [l5l ]. which inter- 
polates between free fermions (N — oo) and the physical 
case of interacting fermions (N = 1): to leading order 
in 1/N, the collision integral with the full medium scat- 
tering T-matrix is consistent with using the density and 
pressure of the free Fermi gas. Specifically, the density 
of a free balanced 2d Fermi gas is 



nX T = 21n(l + z) = 2/9 



(32) 



with thermal length At = \j2it /niksT and fugacity 
z = exp(/3/j) = exp(l/0) — 1 in terms of the reduced 
temperature 9 = T/Tp. The pressure is expressed by the 
polylogarithm Li s (z) as 



P = -nknTOLUl - e 1/e ) 



(33) 



and the internal energy density e — E/V — P equals the 
pressure by scale invariance. The entropy density 



e + P — fin 
f 



= nfc B {-26»Li 2 (l - e 1 ' 8 ) - \n(e l ' e - 1)} 



becomes in the high temperature classical limit 9 — > oo 

s = nk B {2 + ln9 + 0(9- 2 )}. 
Henceforth we will set kg = 1. 

A. Viscosity 

We compute the viscosity of the strongly interacting 
2d Fermi gas with full medium effects. The case with 
Pauli blocking and the bare vacuum scattering cross sec- 
tion, including the limits of high and low temperature, 
has been discussed in Refs. [28l43u| . Our main finding 
is that the medium increases scattering for strong in- 
teraction and thereby substantially lowers the transport 
coefficients, see Fig. Q] For vacuum scattering (squares) 
the system always appears to be in the normal Fermi 
liquid phase and the upturn of the viscosity for low tem- 
peratures is due to Pauli blocking. With medium scat- 
tering the viscosity decreases down to a finite temper- 
ature T c where the medium T-matrix acquires a pole, 
T _1 (g = 0,w = 0) = (Thouless criterion). Below T c 
this pole would formally lead to a diverging collision inte- 
gral C and rj — > in this approximation. A calculation of 
the viscosity in the superfluid B phase of 3 He for T <T C 
found that Pauli blocking and enhanced scattering cancel 
precisely and rj approaches a finite value for T — > . 
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classical gas 
w/ Pauli blocking 
Pauli blocking + medium scattering 
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FIG. 1: (color online) Shear viscosity a — rj/n with and with- 
out medium effects, at strong interaction eb/£f = 2. While 
Pauli blocking (squares) increases the viscosity with respect 
to the classical gas (solid line), medium scattering (circles) 
substantially lowers the minimum as T c is approached from 
above. 



In Fig. [2] the ratio of the viscosity to entropy density 
77/5 is compared for different values of the interaction 
strength. As the binding energy eb is lowered, T c as de- 
fined by the Thouless criterion is shifted to lower temper- 
atures, indicated by the endpoints of the solid lines (the 
endpoints are at T — 1.04 T c ). As an estimate, the min- 
imum for Eb/sf = 0.5 is located at around T/Tp = 0.6 
at a value of rj/s = 0.15, only about twice the proposed 
string theory bound rj/s = l/(47r). 



-52 




FIG. 2: (color online) Viscosity to entropy ratio rj/s with 
medium scattering above T c for different interaction strengths 
eb/ef = 0.1, 0.2, 0.5, 1, 2 (from top to bottom). The 
dashed line indicates the bound l/(47r). 



B. Spin diffusion 

Equivalently, we have carried out the analysis for the 
spin diffusion coefficient D. In the high-temperature 



6 



limit m 



„ QO „ 2 , 2 / 3T 
D = — Q = 7r + In I 



(34) 



the diffusion coefficient depends linearly on 9 with loga- 
rithmic corrections, see Fig. |3l Pauli blocking (squares) 
increases diffusion, while the inclusion of medium effects 
leads to a strong reduction of the diffusion coefficient D 
(circles). 
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3 auli blocking + medium scattering 

1T C ~~ 
T/T F 
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FIG. 3: (color online) Spin diffusion coefficient D in the high- 
temperature limit of a classical gas (solid line) , including Pauli 
blocking (squares) and with the full medium scattering cross 
section (circles). 



V. COMPARISON TO EXPERIMENT 

In order to compare our results for the balanced ho- 
mogeneous 2d Fermi gas with experiments in a trap we 
perform an average over the density profile of the trap 
assuming the local density approximation to hold. At 
high temperatures the density profile in the trap is [28[ 



n(r) 



7 

TT(7 Z 



(35) 



with a 2 = 2T/(mcj^), radial trapping frequency ui± and 
total density J d 2 rn{r) = N. The local Fermi tempera- 
ture is given in terms of the density as 



Tf{t) = —n(r) 
m 

such that the local reduced temperature is 

T mT 



0(r) 



(36) 



(37) 



Tp(r) 7m (r) 
and the local pressure of the free Fermi gas is, cf. (|33|) . 

P(r) = -n(r)T9(r) Li 2 (l - e 1/e(r) ). (38) 



The frequency dependent shear viscosity of the homoge- 
neous system is in kinetic theory jT^j, [H, Eij 



Pt 



1 + L0 2 T 2 



(39) 



in accordance with the viscosity sum rule |45l |. From 
the dimensionless ratio r)(0)/n = a(9) one obtains the 
viscous scattering time 



77(0) 



(40) 



The local viscosity can be defined in terms of the local 
reduced temperature 9(r), 



r)(u,r) 



n(r)a(9(r)) 



l + uj 2 [n{r)a{9(r))/P(r)Y 



(41) 



The spatial integral of the viscosity diverges at u = 
because the dc viscosity is density independent in the 
outer regions of the trap [H, [29[ . In order to obtain a 
finite integral the viscosity is evaluated at the quadrupole 
frequency ujq = \^2oj± [27j |. 



1 

N 



d 2 ri](ujQ,r) . 



(42) 



The global Fermi temperature Tp = v Nuj±_ allows us to 
define a global reduced temperature 8 = T/Tp, so that 
the trap averaged viscosity can be written as 



<«(e» = i / 



d 2 rn(r)- 



a(9(r)) 



1 



'■ a 2 (0(r)) 



(43) 



with dimensionless pressure p(6(r)) = P (r) / (n(r)T) . We 
can change variables and integrate 9(r) = 29 2 . . . oo, 

d0 a{9) 



(a(0)> - 26 2 / 
Ji 



NS 2 p 2 (6) 

Finally, the quadrupole damping rate is [13] 



Q 



2(a(0)) (a(0)) 



NO 



(44) 



(45) 



with (r 2 ) = <t 2 for the density profile in Eq. (|35|) . In 
the high-temperature limit the integrals can be solved 
analytically and yield [28| 



a(9) = 
(a(0)) = 



R9 

27T 7 ' 

m 2 

2tt 

m 



W + In^— j (46) 

r tt 2 a i 
r it 2 n i 

ln [ 1+ ^J (48) 

where we have used p(9) = 1 and (uq/wi) 2 = 2. In 
Fig.Uwe show the quadrupole damping rate vs. interac- 
tion strength and compare with the experimental values 



oj± 2ttVN 
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[27j |. The effect of the medium scattering is most pro- 
nounced at low temperature and strong interaction. This 
leads to strongly enhanced damping, and the peak height 
Tq/w± ~ 0.6 agrees well with experiment, while previous 
theoretical studies found lower peak values Tq/uj± < 0.4 
[29l f30l | . Still, the peak position in our calculation occurs 
at a larger interaction parameter than in the experiment. 
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expt. Vogt et al. (2012) 
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FIG. 4: (color online) Quadrupole damping rate Yq/ui± vs. 
the interaction strength of the trapped gas at T/Tf = 0.3 
and Ef/H = 6.4 kHz, with radial trapping frequency ui± — 
2tt x 125 Hz and N = 2620 particles. 



VI. CONCLUSION 



We have investigated the temperature dependence of 
the shear viscosity and spin diffusion in a two-component 
Fermi gas in two dimensions with contact interactions. 
We used the Boltzmann equation where in contrast to 
former works we took the medium effect due to finite 
fcrmion density into account. We show that the proper 
inclusion of this effect leads to strong suppression of both 
transport quantities. Performing the trap average we find 
that the inclusion of medium effects quantitatively brings 
us rather close to the experimental findings (27| . It is an 
important question for the future to confirm the result 
obtained within the Boltzmann framework with a more 
refined calculation which does not rely on the validity of 
the quasiparticle picture and possibly extends below T c . 
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